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ABSTRACT 

The force free approximation is often useful when describing tenuous plasmas 
in strong cosmic magnetic fields. Time evolution of any such system is gov- 
erned by the information that can be transported along the characteristics of 
the plasma modes allowed by the force free constraint. This brief article eluci- 
dates, for the first time, the nature of the information that can be transported 
along the characteristics for each of the two plasma modes, the Alfven and fast 
modes. It is reassuring that these modes behave as one would expect if the per- 
fect magnetohydrodynamic (MHD) modes in a magnetically dominated plasma 
were evaluated in the limit of zero plasma inertia (i.e., only allowed to propagate 
information consistent with the force free constraint). There are no properties 
of these waves that do not exist within the MHD theory, when this additional 
constraint on the current is imposed. The results of the characteristic analysis 
elucidates the nature of the causality violation in force free black hole magne- 
tospheres. The most significant result from the standpoint of global causality 
is that charge and current perturbations can not be transported along the fast 
mode characteristics in the force free limit. 



1. Introduction 

Applications of force free electrodynamics are common in the astrophysical literature. 
The assumption that the plasma inertia is negligible compared to the magnetic energy density 
is an accurate depiction of the plasma state in many contexts. This condition justifies 
the force free approximation based on a local physical frame analysis. The most notable 
successes of such an approach are descriptions of equilibria in the solar magnetosphere. 
Motivated by the solar work, it has also been extended to studies of the physics of black 
hole magnetospheres (see Blandford and Znajek (1977) and Thorne et al (1986)) and more 
recently gamma ray bursts (see Blandford (2002)). 
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However, the implementation of the force free condition to black hole magnetospheres 
is not so trivial because general relativity is not a local theory. In particular, the concepts 
of an event horizon and ergosphere are globally defined (Hawking 1973). Black hole magne- 
tospheres are the main focus of this treatment, yet all of the results that are derived 
in this article apply equally well to any force free calculation. A salient example 
of the global nature of general relativity is the fact the physical boundary condition at the 
event horizon is that all particles that are near the black hole must corotate with the event 
horizon and go into the black hole on trajectories that approach the ingoing light-like prin- 
cipal null congruence, as seen by any external observer (see Chapter 3 of Punsly (2001)). 
This applies to any plasma, no matter how tenuous and is independent of the magnetic field 
strength. There are no stable equilibria near the horizon, gravity always overwhelms the 
other forces. This property is precisely what distinguishes a black hole from another type 
type of compact object such as a neutron star or white dwarf. All plasma near the horizon 
is inertially dominated: a global condition. In a sense, the plasma inertia is infinite near 
the horizon as viewed globally. By contrast, a local analysis would lead one to believe that 
the force free approximation is valid (i.e., the plasma inertia is neghgible). In spite of this 
defining property of the black hole, force free treatments of black hole magnetospheres that 
extend all the way to the event horizon have been proposed by many authors (most notably 
Blandford and Znajek (1977) and Thorne et al (1986)). 

In this paper, we consider the plasma waves in a force free magnetosphere. This is 
important in order to clarify certain anecdotal defenses of the force free approximation near 
the event horizon which infer that force free plasma waves can transport different information 
than their counterparts in the full perfect magnetohydrodynamic (MHD) treatment of plasma 
waves in magnetically dominated magnetospheres (see Chapter 2 of Punsly (2001) for a 
detailed MHD treatment and Blandford (2001) for the force free discussion). In this article, 
it is shown by means of a straightforward explicit calculation that this claim is erroneous. 
This article elucidates the current structure within the individual wave modes for the first 
time (even the lengthy, complicated and very detailed analysis of certain topics of Uchida 
(1997) in the context of black hole magnetospheres does not address these most relevant 
issues in a force free magnetosphere). 



2. The Force Free Approximation and Magnetically Dominated MHD 

This section compares and contrasts the force free approximation with a magnetically 
dominated perfect MHD system. Generally, the force free approximation is considered as a 
zero mass limit of the full MHD theory. It is usually considered an expedience that reduces 
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the dimension of the full MHD system of hyperbolic differential equations to five. It is very 
rarely argued to be a more accurate physical description than the full MHD theory (but, sec 
Blandford (2002) for such an argument). Generally, the force free approximation is invoked 
is to simplify the mathematical treatment not to improve on the physical depiction. If the 
force free approximation is not pathological then one would expect it to occur as the natural 
limit of zero inertia in the full MHD theory. This section begins to lay the groundwork to 
show that this is indeed the case. The set of force free equations are the current constraint 
written covariantly in terms of the Maxwell tensor and in a local inertial coordinate system, 

F'^V^ = , (2-la) 
PeE + = , (2-lb) 

and Maxwell's equations (written in local coordinates) 

VxE.-if, (2-3) 

V • B = , (2-4) 

V • E = 47rpe ■ (2-5) 

Equation (2.1) also implies the conditions: 

E • B = , (2-6) 
E • J = . (2-7) 

The force free condition in equation (2.1) arises from the zero mass density limit of perfect 
MHD, the plasma has no inertia to stress the magnetic field lines. 



3. Force Free Plasma Waves 

The causal structure of a force free system is governed by the waves supported in the 
force free plasma. In this section we derive the nature of the information that can be 
transported along each of the wave characteristics. Recall that in perfect MHD the Alfvcn 
mode can transport charge and current that was either field aligned (force free) or cross- 
field (inertial). The fast wave cannot transport charge or field aligned currents, only inertial 
currents (see Punsly (2001), Chapter 2). It turns out that the effects of the force free 
constraint are equivalent to the zero mass density limit of perfect MHD plasma waves, as 
expected. Since the fast wave still cannot transport force free current in the zero inertia 
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limit, when the force free constraint on the current is imposed, the only value of current 
consistent with these two concepts (effectively two mathematical equations) is zero current. 
The fast mode in the force free limit carries no current or charge. 

As a consequence of equation (2.6), one can find a frame in which the electric field 
vanishes if the Lorentz invariant F^^'^F^,^ > (i.e., the field is magnetic). There exists a 
coordinate transformation to this frame that resembles a Lorentz boost in the azimuthal 
direction in an axisymmetric magnetosphere except that the boost parameter can exceed 
the speed of light. Thus, this is not a Lorentz transformation in general. However, the angu- 
lar velocity as observed from asymptotic infinity in an asymptotically spherical coordinate 
system is referred to as the field line angular velocity, , even when the azimuthal velocity 
is faster than the speed of light (Blandford 2002; Punsly 2001). Thus, this frame is often 
called the corotating frame of the magnetic field. The hght cyhnder is defined as the sur- 
face at which the azimuthal boost velocity equals the speed of light (in all physical frames) . 
The implementation of this frame facilitates the calculations and one can perform a general 
coordinate transformation in order to evaluate plasma quantities in other physical frames. 
Since, this is not a Lorentz transformation in general, it is best to find coordinate invariant 
expressions. Another frame in which the electric field vanishes is the plasma rest frame. 
However, this frame is poorly defined in the force free limit since there is no accounting of 
the plasma state. Frames in which the electric field vanishes are sometimes called proper 
frames. 

The calculations arc performed in the short wavelength approximation. This is equiva- 
lent to the analysis found in Blandford (2002). No matter what the radius of curvature of 
the magnetic field, one can find plasma mode wavelengths that are short compared to this 
value. This is not very restrictive, even close to the event horizon of a black hole. One should 
be able to find proper frames in which the radius of curvature of the magnetic field is on 
the order of the geometrized mass of a black hole ( 10^^ cm for a supermassive black hole). 
The unperturbed value of the four current, Jq , effectively vanishes in the short wavelength 
approximation that the waves vary on much smaller scales than the background magnetic 
fields (the currents are effectively created by the gradients of the magnetic field). As a re- 
sult of Ampere's law, this approximation is equivalent to treating the magnetic field of the 
unperturbed state as approximately a constant in the proper frame, 

-^B ^ . (3-1) 



We analyze short wavelength linear (small amplitude) perturbations of equations (2.1)- 
(2.7). The waves are chosen in a plane wave representation to have an oscillatory behavior, 
expi(k • r — Lot). We define the x direction to be parallel to k, so all perturbed quantities 
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are Fourier analyzed in {x, t) space. For example, the four current is 

J'^ = + e^(^=^-'^*)(5J'^ (3-2a) 



The interested reader can find a discussion of the effects of field line curvature in a 
black hole magnetosphere in Chapter 6 of Punsly (2001). The nature of the resultant hybrid 
modes is beyond the scope of this treatment that has been chosen to parallel the analysis 
of Blandford (2002). In summary, field line curvature mixes the properties of the plasma 
modes as it does in the earth's magnetosphere (Lanzerotti and Southwood 1979). Long 
wavelength modes are neither fast nor Alfven modes, but have hybrid properties (Appendix 
A is a derivation of the first order effects of magnetic field fine curvature in the force free 
limit) . The faster mode propagates slower as the wavelength increases and is no longer related 
to the critical point structure of the winds that refiect short wavelength mode speeds. The 
implication for black hole magnetospheres is that long wavelength force free fast modes travel 
slower than the speed of light and therefore can not propagate information from near the 
event hroizon to asymptotic infinity. 

The analysis of long wavelength plasma modes is very complicated and one is referred to 
the in depth analysis in Chapter 6 of Punsly (2001) to see the details of outgoing mode prop- 
agation resulting from field line curvature and spacetime curvature near the event horizon 
of a black hole. The results of those calculations are in accord with the causality arguments 
in this article. In particular, one finds that long wavelength plasma modes are extremely 
inefficient in their ability to transport magnetic stresses, field aligned current and charge 
outward from the spacetime near an event horizon. 

In a frame with a vanishing electric field, the perturbed equations (2.1) - (2.7), with the 
condition above become: 

B X 5J = , 

LoSB = -ck X 5B- iAnSJ , 
ujSB — ckx SE , 
k • 5B = , 
ik • (5E = AirSpe , 
5E • B = . 

It is useful to combine (2.2) and (2.3) to get the second order Maxwell's equation, 

V^E - V(V • E) = \ . 



(3-3) 
(3-4) 
(3-5) 
(3-6) 
(3-7) 
(3-8) 



(3-9) 
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This equation gives the following relation amongst perturbed electromagnetic quantities, 

k^SE - k (k • 5E) = + '-^SJ . (3-10) 

c c 

At this point wc explicitly write out the components of the perturbed equation above in a 
local Cartesian basis. The x-axis was previously defined to be parallel to k and we define 
the z-axis such that B always lies in the x-z plane. The components of the second order 
Maxwell's equation become: 

4:7TtiU 

^^y = (^2^2 _ ^2) ' (3-1 la) 

A'jrjijj 
A'jrj. 

SE^ = SJ^ . (3-llc) 



Similarly an expansion of the components of the perturbed force free condition in (3.2) 
yields the following two conditions on the current: 

SJy^O, (3-12) 
SJz — tan 9 SJx , (3-13) 

where 9 is the angle between the propagation vector and the magnetic field direction. There 
is also a third constraint on the current that is obtained by substituting the perturbed second 
order Maxwell's equation (3.11) into the perturbed degeneracy condition (3.8): 

SJjc^k^-uj^)=uj^'^6J,, (3-14) 
^11 

where k\\ and k± are the components of the propagation vector parallel and perpendicular 
to the magnetic field direction, respectively. 

Equation (3.14) can be substituted into the relation (3.13) in order to eliminate 6Jx. 
This yields the dispersion relation for plasma waves in a force free medium, 

[u' - c^e] [uj^ - c^kl] = . (3-15) 



3.1. The Fast Mode 

The dispersion relation (3.15) has two solutions for the phase velocity. The fast mode 
is characterized by the phase velocity, 

^ = ±c . (3-16) 
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The group velocity is 



5B, 


= 0, 


my 


= 0, 




^SEy 




= 0, 


5Ey 




SE, 


= . 



Inserting the value of the phase velocity into the perturbed induction equation, (3.5), and 
the perturbed degeneracy condition, (3.8), along with the vanishing divergence of the mag- 
netic field, (3.6), and the constraint on the perturbed current, (3.14), we can solve for the 
electromagnetic components of the fast wave: 

(3-18a) 
(3-18b) 
(3-18c) 

(3-18d) 
(3-18e) 
(3-18f) 

Substituting the values of the electric field from (3.18d-f) into the second order Maxwell's 
equation, (3.11), and the perturbed Gauss' law, (3.7), the light-like phase velocity of the 
fast wave requires that all of the components of the four-current density that are 
propagated along the characteristics of the force free feist wave vanish identically, 

SJf'^O. (3-19) 

This is a coordinate invariant expression and is therefore true in all coordinate systems. The 
fast mode is linear in the force free limit, since there is a unique wave propagation speed 
given by (3.17). Thus, fast wave packets are linear superpositions of the oscillatory solutions. 
Therefore, the four-current density vanishes for fast wave packets as well. 

Prom (3.18), one also has the following coordinate independent statements: 

^SF^'^SF^^ = , (3-20a) 

= , (3-20b) 

*(5F''"F^, = , (3-20c) 

6F^M^^ = , (3-20d) 

= (3-20e) 

c ° 

where is the group four velocity of the wave, which is expressed in the coordinates of 
the proper frame as: 

. U I) . (3-21) 
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Not all of the relations in (3.20) are independent, for example (3.20d) implies (3.20a). 

The primary result of this analysis, (3.19), that demonstrates the vanishing of the current 
in a fast mode is not a result of a biased interpretation by this author, a miscalculation or 
an unfounded assumption. The same result actually follows from the results of Blandford 
(2002). First of all, the wave structure that is described in a coordinate independent manner 
in (3.20) can be rewritten in the "laboratory frame coordinates" used by Blandford (2002). 
One then finds that, k, 5B and 5E form an orthogonal triad as well as equations (7) and 
(8) of Blandford (2002). These are the results that are used in that treatment to describe 
the causal relevance of a fast wave to a force free magnctosphcrc. Inserting (7) and (8) 
into Ampere's law, either (2.2) or (3.4), yields the fact that 5Jy = and 6Jz = 0. By the 
orthogonality of the triad, SE^ = 0, as a consequence of the definition of the propagation 
direction. Then Ampere's law can be used again to show that 6Jx — 0. Finally, the vanishing 
of the charge density follows from either Gauss' law or current conservation. Thus, (3.19) 
can be found from the analysis of Blandford (2002) quite readily. 



3.2. The Alfven Wave 

The other root of the dispersion relation (3.15) gives the phase velocity of the of the 
Alfven wave: 

^ = ±ccos^. (3-22) 
One can also compute the group velocity, 

du , 'B _ 

Thus, the Alfven mode also propagates information at the speed of light, but only along the 
direction of the magnetic field. 

Combining, the vanishing divergence of the perturbed magnetic field, (3.6), along with 
the constraint on the perturbed current, (3.14), and the induction equation, (3.5), inserted 
into the second order Maxwell's equation, (3.11), one can solve for the components of the 
perturbed electromagnetic field transported along the characteristics of the Alfven mode: 



= , (3-24a) 

ink .4:71 

ckj^ ck± 

SB.^O , (3-24c) 



SBy = -i^SJ, = -i—SJii , (3-24b) 
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5E^ = -i—5J^ = -^^^^11 , (3-24d) 
CO ck ' 

5Ey = Q , (3-24e) 

5E, = - COS eSBy = -i—^5J\\ . (3-24f ) 

ckk±_ 

The component of the current along the magnetic field in the proper frame is defined as 

J|l = J • I . (3-25) 

The perturbed current can be found from (3.22), (3.24) and (3.14), 

cAj I A/ 1 1 II 

5J^ = ^^^^By - -1 J|| , (3-26a) 

SJy^O , (3-26b) 

6J. = ^^6By = -fj, , (3-26C) 

SJu^SJ^i^5By-. (3-26d) 
47r S 

Prom Gauss' law, (3.7), and (3.24) the perturbed charge density is 

5pe = i^SBy = ^ . (3-27) 

The second equality in (3.27) utilized the perturbed law of current conservation which is ob- 
tained from the combination of (3.7) and (3.4). The equations (3.24)-(3.26) can be combined 
to give the following useful relations amongst dynamical quantities that are propagated along 
the Alfven characteristics: 

(5J X B = , (3-28a) 

• B = , (3-28b) 

(5E • 5B = , (3-28c) 

5B • B = . (3-28d) 

Covariantly speaking, (3.28) implies the coordinate independent relations: 



SJ^F^, = , 


(3-29a) 


*6F>'''F^, = , 


(3-29b) 


*6F^''6F^, = , 


(3-29c) 


SF^'-F^, = , 


(3-29d) 


SFi^^SF^, = , 


(3-29e) 


5F^M^^ = , 


(3-29f) 


F,M', = . 


(3-29g) 
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Equations (3.29a-d) are the coordinate independent analogs of the local expressions (3.28a- 
d). From (3.29a), the current flows along the magnetic field lines in a proper frame (this 
expression incorporates the assumption, (3.2), that field line curvature radii are large com- 
pared to the plasma mode wavelength). One can also get the orientation of the propagation 
vector relative to the perturbed and unperturbed fields in all frames by expressing (3.29a) 
in terms of the perturbed covariant version of Ampere's Law, 

SF^'^kaF^, = . (3-30) 

The Alfven mode transports a chcirge and a field aligned current (coveiriantly 
speaking, a force free current). 

4. Causality in Force Free Magnetospheres 

The calculation in the last section allows one to compute the characteristics of short 
wavelength waves in a force free magnetosphere. There are four modes: both ingoing and 
outgoing fast and Alfven waves. The characteristics by themselves are not enough to time 
evolve a force free system, one must also know what information can be propagated along 
these curves. The force free waves carry only two independent pieces of information. In 
order to see this, note that the system has ten parameters: three magnetic field components, 
three electric field components, three current components and the charge. There are eight 
independent relations connecting these quantities - the force free equation (2.1) has three 
constraints, the second order Maxwell's equation, (3.10), has three constraints as well and 
there are the two divergence relations involving the field components in (2.4) and (2.5). This 
leaves two independent quantities that are readily identified from the results of section 3 as 
the parallel (force free) current, Jy and the displacement current orthogonal to the k — B 
plane, SEy. 

4.1. Fast Wave Characteristics 

The fast wave characteristics can be used to propagate finite changes in the displace- 
ment current orthogonal to the k — B plane. Maxwell's equations can be used to find the 
corresponding induced magnetic field. By contrast, no information on the force- free current 

can be propagated along the fast characteristics. If one were to cover a region of spacetime 
with a mesh of short wavelength fast wave characteristics (i.e., the boundary of the spacetime 
region is a essentially a pure fast wave emitter and absorber), and time evolve the system 
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using the method of characteristics then one would find that the current and charge in the 
spacetime would remain unchanged for all time, in all frames. 

There is a claim in Blandford (2001, 2002) that since the fast wave can transport a 
toroidal magnetic perturbation, "it carries information about the poloidal current." This 
conjecture is inconsistent with the nature of the fast wave. The fast wave, like the light wave 
carries only displacement current. The time dependent version of Ampere's law does not 
require a coupling between physical current and the magnetic field. 

4.2. Alfven Wave Characteristics 

By contrast, changes in the parallel current can be propagated along the Alfven wave 
characteristics. However, no changes in in the displacement current orthogonal to the k — B 
plane propagate along these characteristics. From equation (3.27) this current is essentially 
a pure charge moving along the magnetic field lines at the speed of light in the proper frame. 
The Alfven wave can thought of as a charge discontinuity or perturbation propagating at the 
speed of light along the magnetic field (and its associated parallel current) and the resulting 
changes in the electromagnetic field can be found from Maxwell's equations. 

Note that in the degenerate case of wave propagation parallel to the magnetic field, 
there is no k — B plane. In the degenerate case, (3.24) implies that the Alfven wave carries 
no physical current, only a displacement current, 5Ez. 

4-2.1. The Force Free Plasma- Filled Waveguide 

A very important application of the nature of the Alfven wave characteristics is a semi- 
infinite cylindrical waveguide threaded by a uniform axial magnetic field and filled with a 
force free plasma. A rotating conductive disk is suddenly attached to one end of the plasma- 
filled waveguide (see Chapter 2 of Punsly (2001) for pictures and complete details). The 
conductive disk behaves hke a unipolar inductor or Faraday wheel. When the Faraday wheel 
is connected, a Poynting flux begins to be radiated down the plasma-flUed waveguide. The 
Poynting flux is supported by a force free current. In the language of force free waves, the 
Faraday wheel has launched an axisymmetric force free Alfven step wave (see the Appendix 
B for details). 

Note that by (3.26d), the Alfven wave needs to be obhque if it supports a fleld ahgned 
current (i.e., a component of the propagation vector that is orthogonal to the fleld, k±, is 
required). In the plasma- flUed waveguide, the continuity of the tangential electric fleld at 
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the surface of the rotating disk and the degeneracy condition in the force free plasma imply 
that there is a poloidal electric field, E^, that is launched by the Faraday wheel orthogonal 
to the poloidal magnetic field, , 

= . (4-1) 

c 

The coordinate, p, is the radial coordinate in a cylindrical coordinate system and the field line 
angular velocity, Q^. is identical to the disk angular velocity, . The radial variation of the 
electric field in (4.1) and in the oscillatory modes of (Bl) is tantamount to a standing wave 
component of the propagation vector that is orthogonal to the magnetic field (note that the 
same wave polarization phenomenon occurs in vacuum electromagnetic waveguides). The 
step Alfven wave in the cylinder (constructed in Appendix B) is a superposition (Fourier 
integral) of the oblique, axisymmetric Alfven normal modes of the waveguide. The standing 
wave properties (obliquity) of the axisymmetric normal modes in the cylindrical waveguide 
are a consequence of the boundary conditions on the cylinder walls. 



4-2.2. Pulsar Magnetospheres 

One would expect a strong analogy between the plasma-filled waveguide and an ideal- 
ized, force free, field aligned pulsar magnetosphcrc. The neutron star is the unipolar inductor, 
the magnetosphere is the plasma-filled waveguide and the wind is the Poynting flux in this 
analogy. Thus, one would expect that the pulsar wind is essentially an Alfven wave radiated 
by the neutron star in this idealized limit. 

Unfortunately, the causal nature of MHD in a relativistic magnetosphere was obfuscated 
by anecdotal comments in Blandford (2001) through an analysis in the force free approxima- 
tion. If the force free analysis is not pathological then it must agree with the magnetically 
dominated limit of MHD. There is a claim that for axisymmetric modes, the Alfven wave 
can only support variations in the poloidal magnetic field, not toroidal variations. However, 
this cannot be correct as the example of the plasma-filled waveguide shows by the explicit 
construction of the axisymmetric Alfven modes that transport toroidal magnetic fields. The 
concept that a standing wave component of the propagation vector that is orthogonal to the 
magnetic field can exist in an axisymmetric mode (as in the plasma-filled waveguide) was 
not understood in Blandford (2001). The existence of such a wave vector polarization im- 
plies, with the aid of (3.30) that axisymmetric Alfven modes that are capable of transporting 
toroidal field perturbations along the direction of poloidal magnetic field most certainly exist 
in a pulsar magnetosphere. This oversight allows for the claim that the force free analysis 
leads to an "opposite conclusion" to that of the MHD analysis. In reahty, if one is prudent. 
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it can be shown as it has be done here that the two analyses agree, as they must. 

5. Applications to Black Hole Magnetospheres 

The force free approximation has been used in models of black hole magnetospheres in 
which spin energy is extracted electromagnetically from the black hole (Blandford and Znajek 
1977; Thorne et al 1986). The force free wind systems in these treatments are characterized 
by two constants in each magnetic flux tube: the poloidal current (equivalently the toroidal 
magnetic field or angular momentum flux) and the field line angular velocity (this is quantity 
is equivalent to specifying the ratio of energy flux to angular momentum flux, the Goldreich- 
Julian charge density and the cross field potential). 

The mathematical method of solution implements the force free condition at the horizon 
as a boundary condition that is used to determine the field aligned current density and the 
field line angular velocity in each flux tube. Since, by (3.17), the fast wave can propagate in 
any direction, at the speed of light, outgoing fast waves can propagate from anywhere outside 
of the event horizon (the fast critical surface for outgoing waves is the event horizon). By 
contrast, the "light cylinder" located within the ergosphere is the Alfven critical surface, 
which no outgoing Alfven wave can traverse (Blandford 2002; Punsly 2001). For rapidly 
rotating black holes, this surface is sufficiently far from the event horizon that the constraints 
imposed by the spacetime metric near the horizon are irrelevant. Thus, the use of the 
horizon boundary condition can only be justified from a causal perspective if the fast mode 
is involved in determining the wind constants. However, by (3.19), changes in the poloidal 
current or charge can not be transmitted along the fast wave characteristics in the black hole 
magnctosphcre. Thus, the spacetime near event horizon can not affect changes in the wind 
parameters. The use of the event horizon as a causal wind boundary is not permitted in 
the force free approximation. Consequently, the mathematical method of solution is flawed. 
The same causal structure is known to exist in magnetically dominated perfect MHD (see 
Chapter 6 of Punsly (2001)). Going to the force free hmit was conjectured in (Blandford 
2001, 2002) to obviate any need for concern with this method of solution that were revealed 
in the full MHD analysis. This claim is clearly not supported by the results of this article. 

The causal structure of the black hole magnetosphere in the force free limit has been 
shown, in this article to behave formally as the zero mass limit of MHD. As in perfect MHD, 
the wind parameters must be created outside of (farther out from the black hole than) the 
ergospheric Alfven critical surface. 

A further critique of the MHD analysis of black hole causahty was that the Alfven 
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wave really does not transport changes in current and charge, only oscillations (Blandford 
2001). However, this statement is not an accurate depiction of the nature of propagating 
disturbances in MHD. It is well known that small amplitude MHD step waves obey the 
same dispersion relations as oscillatory waves. Alfven step waves and Alfven wave packets 
in general obey the same constitutive relations amongst field and dynamical quantities and 
propagate along the same characteristics as the oscillatory waves (Kantrowitz and Petschek 
1966). This also follows by taking the appropriate Fourier integral to create a step wave 
from sinusoidal waves. The Alfven wave packets can transport finite changes in the charge 
density and field aligned currents in a force free magnetosphere. This is shown by explicit 
construction in Appendix B. 



6. Discussion 

There are three main reasons for not considering black hole magnetospheres that are 
force free everywhere to be physically realizable: 

1. As mentioned in the Introduction, the physical plasma state near the event horizon 
is always inertially dominated. This physical state is exactly the opposite of the zero 
inertia limit that is inherent to the force free assumption. 

2. Prom equation (2.7) there is no region within the force free magnetosphere where 
energy fiux (Poynting fiux) and angular momentum fiux can be created. Energy and 
angular momentum fiux (essentially, the poloidal current and Goldreich- Julian charge 
density) must be injected into the wind system at the boundaries. If the magnetosphere 
is force free everywhere then the energy flux (current and charge) must emerge from 
the horizon. However, the charge and current necessary to support the energy and 
angular momentum flux can not be propagated along the plasma wave characteristics 
emanating from the spacetime near the horizon as discussed in the last section. Thus, 
such a wind system is inherently acausal. 

3. There are three MHD computer simulations of magnetically dominated black hole 
magnetospheres (the magnetic energy density exceeds the plasma energy density in the 
initial state) that have been published. All of them evolve to a state with strong cross- 
field (incrtial, not force free) current fiow in the ergosphere (Camenzind and Khanna 
2000; Koide et al 2002; Semenov et al 2001). The inertial current fiow initiates in 
regions far enough from the event horizon that it is not interpretable as a horizon 
surface current. The strong cross- field current is a necessary component for extracting 
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black hole energy in the theory of black hole gravitohydromagnetics (GHM) as well in 
Punsly (2001). 

It should be noted that a wind system beyond (farther from the black hole than) the 
ergospheric dynamo region for the toroidal magnetic field (poloidal current) could in prin- 
ciple be well described by the force free approximation. From the nature of the force free 
characteristics, this dynamo region can exist only in the spacetime farther from the black 
hole than the ergospheric hght cylinder (Alfven critical surface). The Poynting flux can 
propagate outward from the dynamo region to asymptotic infinity as an essentially force free 
wind. Such a wind system is described in detail in Chapter 9 of Punsly (2001). 

I would like to thank Ferd Coroniti for his valuable and intelligent insights. I am also 
indebted to Serguei Komissarov for his useful comments regarding force free plasmas and 
propagating discontiuities. 



A. First Order Curvature Effects - Hybrid Modes 

In this appendix, we calculate the first order effects of field fine curvature on the force 
free small amplitude plasma modes in a proper frame. One finds new hybrid modes as in 
MHD. The first order corrections to the short wavelength approximation can be found in 
the eikonal approximation. The sinusoidal solutions described in (3.2) are replaced by: 

where k is a slowly varying vector valued function. In the present circumstance, the unper- 
turbed current no longer vanishes by assumption. However, by the force free condition, it 
flows parallel to the magnetic fleld in the proper frame. This can be expressed in terms of a 
slowly varying "radius of curvature" , i?, as 

J„ = ci?-^B„ . (A2) 

The validity of the eikonal approximation in (Al) requires that: 

I k I i? > 1 , (A3) 



The perturbed Maxwell's equations and the perturbed degeneracy condition carryover 
from the sinusoidal case, (3.4) - (3.11), unchanged. The only perturbation equation that 
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changes is the force free condition, (3.2), that now becomes: 

X 5J + X Jo = . (A5) 

The perturbed force free relation found from expanding the "x" and "z" components of (A5) 
is 

^^^.{kR^SJ^^SJ,. (A6) 

Note that by (A3), the large radius of curvature approximation, equation (A6) is approxi- 
mately (3.12). The "y" component of the perturbed force free relation (A5) combined with 
(A6) yields the relation: 

[C2A;2_^2]2 k " ' k ' ^ ^ 

Note that in the hmit of a large magnetic field radius of curvature that (A7) becomes (3.13). 
Thus, the complete set of perturbed equations reduce to the equations found in section 3 as 
the curvature of the field gets small. 

The dispersion relation is found by combining (A7) with the perturbed degeneracy 
condition (3.14), 

y - c'^e] y - c^A^J] = 16n^eklc\kR)-^ . (A8) 
There are two roots of the dispersion relation.: 



c^k^ + \^ [^A;j + 647r2A;2;t2(A;i?)-2 - k\^ , (A9) 
A;l + 647r2HF(A;i?)-2- fci] . (AlO) 



a;2 _ c^itj - |c2 

Differentiating (A9) and (AlO) yields the group velocities of the two modes. For parallel 
propagation, the group velocity is the speed of light up to terms that are O [(/ci?)"^]. The 
parallel propagation assumption is invalidated at second order since the magnetic field is 
curved. In the limit that [k.^^Kf' ^ 1, one can obtain the following Taylor series expansions 
of the group velocity for oblique propagation: 

'^'^ ■ ±c[l-87r2(A;i?)-2cot2^] , (All) 



^ 0+167r2cot2^(A;i?)-2 

dcu _ , ccos^ 

\/l-167r2 cot^ e{kR)-^ 



^^±— ^==^^====^±0 cose [l + 87r\kR)-^ esc'' e] . (A12) 



Equations (A9) and (All) approximate to (3.16) and (3.17) as the radius of curvature 
becomes infinite. Thus, these equations represent the fast hybrid mode. Similarly, (AlO) 
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and (A12) become (3.22) and (3.23) in the short wavelength hmit. These equations represent 
the hybrid Alfven wave. Notice that the fast hybrid mode propagates slower than the speed 
of light for oblique waves propagating across a curved background magnetic field. From 
(A12), the Alfven mode does not propagate perpendicular to the field. Thus, there is not 
Alfvenic properties for the fast mode to couple into this case. Thus, the hybrid fast mode 
still propagates at the speed of light for perpendicular propagation as indicated in (All). 

There are implications of this result to the primary focus of this article. Since the 
fast wave no longer propagates isotropically, it can now carry a physical current (not just 
a displacement current). Denote the displacement current as, J^. Then (3.11a), (A6) and 
(A9) imply that the ratio of displacement current to physical current for oblique fast hybrid 
waves is: 



(A13) 

J|| 647r ^ ' 

Consequently, the ratio of displacement current to the physical current diverges in the short 
wavelength limit as in (3.19). Similarly, from the perturbed Gauss' law (3.7a), (3.11c) and 
(A9), one can find the ratio of displacement current to the charge density: 

— ~ c- — . (A14) 

Pe 647r ^ ' 

Comparing (A9) to (A13) and (A14), note that as the fast wave propagation speed slows 
down due to field line curvature, it carries more physical current and charge. It becomes 
increasingly more like the Alfven wave as the field line curvature increases. However, (A13) 
and (A14) show that in the eikonal approximation, the fast wave is still almost purely 
electrodynamic: a propagating displacement current with a minute physical current that is 
only O [{kR)~^] of the displacement current. 



B. The Alfven Wave Packet in the CyUndrical Waveguide 

In this appendix it is shown by explicit construction that a superposition of oscillatory 
Alfven waves (an Alfven wave packet) does indeed transport field aligned currents and charge, 
not just oscillatory perturbations. We rely on the magnetically dominated MHD treatment 
of the plasma-filled cylindrical waveguide found in Chapter 2 of Punsly (2001) and pass to 
the force-free limit. In order to save space, we do not reproduce the results of Chapter 2 here, 
but refer to them in detail and the reader is urged to consult that reference. In summary, 
a cylindrical waveguide filled with a tenuous plasma is placed inside of an infinitely long 
solenoid so that it is threaded by a strong axial magnetic field. The waveguide is semi- 
infinite in the sense that one end is terminated by a rotating conductive wheel. This rotating 
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conductor or Faraday wheel (also known as a homopolar generator or unipolar inductor) is 
an active element in the waveguide circuit that radiates plasma waves. In the first subsection, 
we construct the relevant axisymmetric, m=0, perfect MHD Alfven modes in the cylinder 
using the detailed results of Punsly (2001). In the next subsection, we compare and contrast 
a wave packet (a step wave) created from a linear superposition of the previously calculated 
oscillatory Alfven modes to the radiation emanating from the Faraday wheel (the nature 
of this radiation was derived previously in terms of a formal MHD wind solution in the 
waveguide in Punsly (2001)). Finally, it is noted that the depiction of the wind as a pure 
Alfven step wave becomes increasingly more accurate as one passes to the force- free limit. 



B.l. Components of the Oscillatory Waves 

This section constructs an oscillatory solution from Chapter 2 of Punsly (2001) for an 
axisymmetric Alfven wave, m = 0, in the plasma-filled waveguide that is consistent with the 
attachment of the Faraday wheel (in the realm of magnetically dominated perfect MHD). 



B.1.1. The Electric Field 

From equation (2.93b) of Punsly (2001), the radial electric field in an axisymmetric 
Alfven wave in a cylindrical waveguide is arbitrary and is determined by the boundaries (the 
field lines can shear relative to each other, hence these modes are often called "shear Alfven 
waves"). Considering the frozen-in condition within the Faraday wheel and the continuity 
of the tangential electric field at the Faraday wheel/plasma interface let us choose the radial 
electric field in the Alfven wave to be (in cyhndrical coordinates): 

The other electric field components are given in Punsly (2001) by (2.86c) and (2.93a) as 

(B2) 

E' ^0. (B3) 



B.1.2. The Currents 
From equations (B1)-(B3) and Gauss' law, the charge density is; 



Pe = Jh&Sk.---^t) _ /B4N 

27rc 
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Prom the Alfven wave dispersion relation, (2.82) of Punsly (2001), and the frozen-in form of 
the momentum equation, (2.86a) of Punsly (2001), along with equation(Al), the cross-field 
current density is 

jp ^ ^ e^(fe--t) , (B5) 

A 

where is the intermediate three speed and t/^ is the pure Alfven speed. In the force free 
limit, the pure Alfven speed diverges as the intermediate speed approaches the speed of light. 
Prom the law of current conservation, 

^ + V-J = 0, (B6) 
combined with (B4) and (B5) yields the axial current density: 



B.1.3. The Magnetic Field 

Using the value of the axial current in (B7) in Ampere's law yields the toroidal magnetic 
field strength: 

The other components are given by (2.92a) and (2.96) of Punsly (2001): 

5^ = , (B9) 
B" . (BIO) 



B.2. Making Wave Packets 

We can construct wave packets of the oscillatory solutions above by taking Pourier 
integrals. We are interested in the step wave. Por the Alfven wave, one can use the dispersion 
relation (2.82) of Punsly (2001) in the cyhnder to write, 

gj(fez2-a;t) _ ^-ikzivjt-z) (Bll) 

The step function is given by. 



e{vjt -z) = 




(B12) 
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We construct a wave packet of Alfven waves, with a spectral amplitude, A{kz), 

/ + 00 
AiQiPiQdk, , (B13) 
-oo 

where ip{kz) is one of the oscillatory wave function components in equations (Bl)-(BIO) and 
e is an arbitrarily small positive number in the usual sense. The wave packet of oscillatory 
solutions with field and current components given by (Bl)-(BIO) and spectral amplitude 
given by (B14) has the following field and current distributions (for the sake of demonstrating 
equivalence, the corresponding equation number from Punsly (2001) that were derived by 
the MHD wind calculation in the cylinder is placed on the right hand side of the equals sign 
in these relations): 



E" = !^B^Q{v.t -z)^ (2.98) , (B15) 

c 

= , (B16) 

E' , (B17) 

Pe = -^^Q{v,t -z) = (2.99) , (B18) 

A 

r = ^c" = -^£^e{vA -z)^ (2.119) , (B20) 

Vj 2'KVj 

B<l> ^ _^SEpQ(y t - z) = (2.118) , (B21) 

= , (B22) 

BP^O. (B23) 



These quantities agree with the downstream state of the plasma that was found near the 
flow front in Punsly (2001) in the limit that — Vj. Notice that Jf vanishes downstream 
of the wavefront and only has a surface component on the wavefront as it does in equation 
(2.110) of Punsly (2001). 



B.3. Conclusion 



Prom equations (BIS) and (B20), the Alfven wave packet actually transports charge 
and fleld aligned current. 
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The solution is not exact. There are errors associated with the small inertial terms in 
the magnetically dominated limit. At the wavefront, the relativistic MHD shock equations 
do not solve exactly. These equations are the conservation of the axial components of the 
stress energy tensor across the discontinuity in the frame of the wavefront. Pressure balance 
(magnetic pressure from the toroidal magnetic field in the frame of the wavefront) can not be 
achieved because the Alfven wave, unlike the fast wave, has no compressive properties. The 
errors in the stress-energy balance can be found in the frame of the propagating discontinuity 
by using (B15) and (B21) to compute the proper toroidal magnetic field. The errors in the 
shock equations are on the order of (Jl^pC/~^)^. This is biquadratic in two small quantities: 
fi^p/c <S 1 by construction and by the magnetically dominated condition, cU~^ <C 1. 
The violation of the shock relations are therefore a very small second order effect in the 
magnetically dominated limit. Notice that the errors vanish completely in the force-free limit 
as the Alfven three speed approaches the speed of light. In magnetically dominated perfect 
MHD, these small errors are accounted for by a fast switch-on shock that creates a and 
an E"^ downstream of the shock front. This is an MHD precursor, an infinitesimal distance 
upstream, to the Alfven rotational discontinuity at the terminus of the Alfven step wave. As 
demonstrated by explicit construction above, it is the Alfven rotational discontinuity that 
imprints the charge and field aligned current on the waveguide plasma. The interpretation 
of the waveguide wind solution as an Alfven wave is exact in the forcc-frcc limit and is 
extremely accurate to first order in the magnetically dominated MHD limit. However, in 
relativistic MHD it is important to note that both fast and Alfven modes are needed to solve 
this type of "Riemann problem." 

It was an oversight in Punsly (2001) that nonlinear fast waves such as an abrupt dis- 
continuity can in principle transport charge and current. However, in practice, due to the 
limited polarization properties of such waves (the fast discontinuity can only affect changes 
in the magnetic field in the plane spanned by the normal to the wavefront and the magnetic 
field upstream) unless the magnetic field is perfectly homogenous (not the case in a magneto- 
sphere) they must be accompanied by an Alfven rotational discontinuity in a magnetosphere 
in order for the discontinuities to be solutions of the relativistic shock equations (Kantrowitz 
and Petschek 1966). The important point is that in general both modes are needed. The 
discussions of Punsly (2001) were designed to show that the Alfven wave is always required 
in MHD, so that the fast wave cannot determine the wind constants alone - and this conclu- 
sion is unchanged and hence the MHD causality arguments are unchanged. In the force-free 
limit, the situation is more extreme. Since the fast wave is linear in the force free limit (the 
fast speed is always c), all wave packets constructed from the complete set of oscillatory fast 
modes will also carry no current or charge in any frame (see section 3.1). Thus, the Alfven 
wave is responsible for setting the wind constants in a force free black hole magnetosphere. 
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This discussion corroborates the physical accuracy of the assumptions used in the one 
dimensional string analysis of flux tubes in Semenov et al (2001). They actually ignore the 
fast mode completely (because as they say, the fast wave is well known to primarily provide 
only pressure changes in an MHD fiux tube) and only consider the Alfven mode. This seems 
like a very reasonable assumption based on the analysis provided here. 
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